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Abstract 

Ac susceptibility measurements were performed on discontinu- 
ous magnetic multilayers [Co8oFe2o(t)/Al203(3nm)]io, t = 0.9 and 
1.0 nm, by Superconducting Quantum Interference Device (SQUID) 
magnetometry. The CoFe forms nearly spherical ferromagnetic single- 
domain nanoparticles in the diamagnetic AI2O3 matrix. Due to dipo- 
lar interactions and random distribution of anisotropy axes the system 
exhibits a spin-glass phase. We measured the ac susceptibility as a 
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function of temperature 20 < T < 100 K at different dc fields and as 
a function of frequency 0.01 < / < 1000 Hz. The spectral data were 
successfully analysed by use of the phenomenological Cole-Cole model, 
giving a power-law temperature dependence of the characteristic re- 
laxation time Tc and a high value for the polydispersivity exponent, 
a f« 0.8, typical of spin glass systems. 
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INTRODUCTION 

The dynamic and static magnetic properties of spin glasses (SG) are still 
a subject of intense experimental and theoretical research. In the field of 
experiments a vast variety of different spin glass systems ( [Young, 1997] ) have 
yet been found and investigated. A rather new class are the so-called su- 
per spin glass (SSG) systems ( Kleemann, Fetracic, Binek, Kakazei et al. 



200 1| ). Here the sample is composed of an ensemble of dipolarly interact- 



ing nanoparticles, each having a superspin moment in the order of 1000 fiB- 
The spin glass properties are due to frustration, a natural property of dipolar 
interaction, and to randomness of the anisotropy axis directions, frequently 
also of the spin sizes. Two different types of realisations of SSG systems 
exist, frozen ferrofiuids ( pormann, Fiorani and Tronc, 19"9^ ; [Dormann, Fio- 



rani, Cherkaoui, Tronc et al, 1999| ; PJurberg, Svedlindh, Nordblad, Hansen 



et al, 1997| ; [Mamiya, Nakatani and Furubayashi, 199"9|) and discontinuous 
magnetic multilayers ( [Sankar, Pender, Borchers, Smith et al., 200U| ; [Sousa 



Kakazei, Pogorelov, Santos et al., 2001 ; Kleemann, Petracic, Binek, Kakazei 



et al., 200T| [Petracic, Kleemann, Binek, Kakazei et al., 2002| ; [Sahoo, Petracic 



Binek, Kleemann et al, 2002a[ ). 

It is widely accepted that 3-dimensional (3D) nanoparticle systems with 
high enough density of the particles and sufficiently narrow particle size dis- 
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tribution do have spin glass properties, i.e. there exists a phase transition 
temperature, Tg, where the characteristic relaxation time and the static non- 
linear susceptibility diverge ( pjurberg, Svedlindh, JNordblad, Hansen et al. 



1997] ; Kleemann, Petracic, Binek, Kakazei et al, 2001 ; Petracic, Kleemann 



Binek, Kakazei et al, 200^ ; |Sahoo, Petracic, Binek, Kleemann et ai, 2002aj) . 



In order to observe SSG properties the collective glass temperature, Tg, has to 
be larger than the so-called blocking temperature, Tf,, at which the relaxation 
time of the individual moments (|Neel, 19491 ; [Brown, 19631) 



T = Toexp{KV/kBT), (1) 

reaches the order of the timescale of the experiment. Here K is the effective 
anisotropy constant of one nanoparticle, V the volume and tq ~ lO^^^s the 
relaxation time at T ^ oo. Below the particle moments are "blocked". 

The condition Tg < T^ is met in our discontinuous metal-insulator mul- 
tilayers (DMIMs) [Co8oFe2o(t)/Al203(3nm)]„, where t < 1.0 nm is the nom- 
inal thickness of the ferromagnetic CoFe layers and n the number of bilay- 
ers. The CoFe does not form a continuous layer but forms nearly spherical 
particles embedded in the diamagnetic AI2O3 matrix ( [Kakazei, PogorelovJ 
Lopes, Sousa et al, 200"T[) . One finds self-organized arrangements of parti- 



cles in each layer ([Stappert, Dumpich, Sahoo, Petracic et al., '200'2 ), i.e. the 



inter-particle distances are nearly constant. While SSG behavior is found for 
relatively small values of t < 1.0 nm and n = 10 ( [Kleemann, Petracic, Bine¥ 



Kakazei et al., 200 1[; [Petracic, Kleemann, Binek, Kakazei et al., 2002| [Sahoo 



Petracic, Binek, Kleemann et al, 2002a[) , for higher values of the nominal 
thickness, 1.0 < t < 1.4 nm, superferromagnetism (SFM) is observed ( [Klee- 



mann, Petracic, Binek, Kakazei et ai, 20"0T| [Sahoo, Sichelschmidt, Petracic 



Binek et ai, 2002b| ; [Chen, Sichelschmidt, Kleemann, Petracic et ai, 2002[) . 



In this article we will focus on the SSG systems [Co8oFe2o(i)/ 
Al203(3nm)]io, with t = 0.9 and 1.0 nm. The existence of a spin glass 
phase was evidenced by means of dynamic criticality, static criticality of the 
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non-linear susceptibililty and dynamical scaling (|Kleemann, Petracic, Binek 



Kakazei et ai, 200 1| ; Petracic, Kleemann, Binek, Kakazei et al, 2002| [Salioo 



Petracic, Binek, Kleemann et al., 2002a| ). All three methods yield convincing 
values for the glass transition temperature and the dynamical critical expo- 
nents, respectively, Tg 44 K, zu ^ 9.5, 7 ~ 1.47 and (3 ~ 1.0 for t = 0.9 nm 
and Tg ^ 49 K, zu ^ 10.0, 7 ^ 1.36 and f3 ^ 0.6 for t = 1.0 nm. Tg and zu 
are the error weighted average values obtained from different methods. 

While zu characterizes the divergence of the relaxation time of the largest 
ordered cluster as T Tg, there is a wide distribution of shorter relaxation 
times due to non-percolating clusters. They are characteristic of the glassy 
nature of the system and deserve a focused investigation, which will be de- 
scribed in the present paper. To this end we analyse the results of measure- 
ments of the complex ac susceptibility carried out at different ac amplitudes 
and bias fields and frequencies /. In particular the so-called Cole-Cole pre- 
sentation, x" vs x'5 will be discussed in terms of appropriate empirical models 
of relaxational polydispersivity. 



EXPERIMENTAL 

The DMIM samples Glass/Al203(3nm)/[Co8oFe2o(t)/ Al203(3nm)]io {t = 
0.9 and 1.0 nm) are prepared by sequential Xe ion beam sputtering from two 
seperate targets ([Kakazei, Pogorelov, Lopes, Sousa et ai, 200 1|) . The CoFe 
forms nearly spherical granules of approximately 3 nm diameter and 2 nm 
inter-particle spacing as found from transmission electron microscopy (TEM) 
studies ( ptappert, Dumpich, Sahoo, Petracic et ai, 2002|) . 

The measurements were performed by use of a commercial Supercon- 
ducting Quantum Interference Device (SQUID) magnetometer (MPMS-5S, 
Quantum Design). The ac susceptibility, X = x' ~ ^x" ^ is extracted from 
the linear response of the sample on an oscilUating ac field, /io-f^ac = 0.05 or 
0.4 mT at different ac frequencies, 0.01 < / < 1000 Hz. The constant dc 
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field was either fioH = (0 ± 0.03) mT or (0.6 ± 0.1) mT. 



RESULTS AND DISCUSSION 

Figure |l| shows the real x' and the imaginary parts x" of the ac susceptibility 
vs temperature T for the samples t = 0.9 nm (a) and 1.0 nm (b) under 
four different conditions. Curves 1 and 1' are measured at the ac frequency 
/ = 0.1 and curves 2 and 2' at 1 Hz, whereas for curves 1 and 2 an ac field 
amplitude of fioHac = 0.05 mT and a dc field of fioH = mT were applied. 
For curves 1' and 2' an ac field amplitude of /io-f^ac = 0.4 mT and a dc field 
of fioH = 0.6 mT were used (see Figure ^ for an illustration). For both 
samples a similar behavior is encountered. Both the increase of the probing 
ac field amplitude and the application of a bias field result in a supression of 
the amplitude of the real part x'{T) and a shift AT^ of the peak to higher 
temperatures. Quantitatively the shift is ATm(l' — 1) = 2.2 (3.8) K and 
ATm(2' — 2) = 2.3 (3.6) K for t = 0.9 (1.0) nm, respectively. The imaginary 
part x"(T) is also suppressed, but the infiection point at Tj is shifted to 
lower temperatures, Ar/(1' - 1) = -6.0 (-10.1) K and AT^(2' - 2) = -6.2 
(—8.2) K for t = 0.9 (1.0) nm, respectively. This behavior is well known 
from other SG systems and model calculations ( Panella and Mydosh, 19"?^ ; 



Barbara, Malozemoff" and Imry, 1981D and can be explained in terms of a 



competition between the non-critical linear susceptibility and the critical 
non-linear susceptibility. In other words, the suppression of both the real and 
the imaginary parts refiects the obvious fact that the M{H) curve becomes 
increasingly non-linear when increasing the ac amplitude and/or the bias 
field (see Fig. |). 

Next we studied the frequency spectra, x'if) x"{f)- Figure ^ and 
Fig. I (t = 0.9 and 1.0 nm respectively) show the real x' (a) aiid the imaginary 
part x" (b) as functions of the ac frequency / for different temperatures T = 
45, 50, 55 and 60 K in zero-field and HoHac = 0.05 mT. While some negative 
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curvature still indicates a well-defined dispersion step at / > 10^ Hz for 
T > 60 K, this step becomes gradually broadened as T decreases. At low T 
the real parts show nearly constant negative slopes, thus corresponding to an 
extremely broad dispersion step. The imaginary parts reveal extremely broad 
peaks, which strongly shift to lower frequencies with decreasing temperature. 
Obviously our SSG system exhibits a very wide distribution of relaxation 
times with a pronounced temperature dependence. 

A satisfactory description of the data is provided by the phenomenolog- 
ical Cole-Cole model ( |Cole and Cole, 1941| ; [Jonscher, 1983| ) and was suc- 
cessfully applied e.g. to 2-dimensional (2D) ( pekker. Arts, de Wijn and 



van Duyneveldt, 1989| ; Pagiwara, 1998| ) or pseudo-l-dimensional SG systems 



( Ravindran, Rubenacker, Haines and Drumheller, 1989|) . The complex ac 



susceptibility, X = x' ~ ix" i is written in the Cole-Cole model as ([Jonscher ,| 
198^ ) 



where xo and xs are the isothermal (low-f) and adiabatic (high-f) suscepti- 
bilities, respectively, t^. is the characteristic relaxation time and a a measure 
of the polydispersivity of the system. The case a = yields the standard 
Debye-type relaxator with one single relaxation frequency, as found, e.g., in 
the case of a monodisperse ensemble of non-interacting superparamagnetic 
particles obeying Eq. |T[ The limiting case a = 1 corresponds to an infinitely 
wide distribution of relaxation times. In SG systems one expects values of a 
near to 1. 

After decomposing Eq. ^ into its real and imaginary parts it is possible 
to perform a fit to the data as shown in Fig. ^ and ^. One finds (compare to 
Dekker et ai, 1989; Ravindran et al, 1989) 



// N , Xo-Xs (-. sinh[(l - a) In(ojrc)] \ 

X ^ = Xs + — ^ — 1 7771 vTT ^TT ^71 ^ 

2 \ cosh[(l — «) ln(c<jrc)J + cos[^(l — «)7rjy 
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2 Vcosh[(l-a)ln(cure)] + cos[i(l-a)7r]y ' 
where u = 2iTf. Best results are obtained, when fitting to the imaginary 

part x"{f)^ since only three parameters, Xo ^ Xs, (1 ~ ct) and Tc are needed 

in this case. 

Figure | shows the results from the fitting, Tc (open circles) and a vs T 
(open diamonds) for both samples, t = 0.9 (a) and 1.0 nm (b). One finds that 
the characteristic relaxation time Tc is increasing with decreasing tempera- 
ture. It changes by eight (a) or ten orders (b) of magnitude, respectively. By 
this kind of extraction of Tc one has access to an extremely wide timescale and 
is, hence, more advantageous compared to the standard method of extracting 
Tc from the x'{T) data. It is straightforward to perform a fit of the Tc(T) 
data to a critical power-law, which was already used in previous publica- 
tions ([Kleemann, Petracic, Binek, Kakazei et ai, 2001| ; Petracic, Kleemann 



Binek, Kakazei et al., 200^ ), Tc = tq{T /Tg — 1) ^'^ (solid line). It yields rea- 
sonable results, but the value of zv must be kept restricted or even fixed to 
zu = Q. Then we obtain To = (5.0 ■ lO^^ ± 5.1 • 10-^) s,Tg = (42.63 ±0.18) K 
and zv = 9.0 ± 0.7 (t = 0.9 nm) and tq = (3.46 • 10"^ ± 2.1 • 10^^°) s, 
Tg = (43.9354 ± 0.0002) K and zu = 9 fixed {t = 1.0 nm), respectively. 
In the case of t = 0.9 nm the values for Tg and tq correspond well to the 
values obtained previously ( [Kleemann, Petracic, Binek, Kakazei et al., 2001 



Petracic, Kleemann, Binek, Kakazei et at., 2002|) . This does not apply to the 



t = 1.0 nm sample, where Tg ^ AA K differs strongly from the value shown 
above, Tg ^ 49 K. Interestingly the Cole-Cole fit to the T = 45 K data for 
t = 1.0 nm does not converge (encircled data points in Fig. |^ (b)) leading to 
the conclusion that the data emerge from the non-ergodic regime, T < Tg. 
It is worth to mention that the fit to the modified power law according to 
QSouletie and Tholence, 1985| ), = tq{1 — Tg/T)~^'-', (broken line) yields 



similar values, i.e. tq = (2.67 • 10"^ ± 2.9 • 10"^) s, Tg = (43.19 ± 0.11) K 
and zu = 9.0 fixed {t = 0.9 nm) and tq = (3.4 ■ 10"^° ± 6.2 ■ 10"^°) s, 
Tg = (44.004 ± 0.036) K and zu = 10 fixed (t = 1.0 nm), respectively. It is 



not possible to judge about the advantage of this method here. 

The exponent a increases, as expected, with decreasing temperature 
(Fig. ||. Its high value [a ~ 0.8) meets the expectation, that a SG system 
should have a very broad distribution of relaxation times QMydosh, 19931) . 

Often suscpetibility data are presented in a way, where the imaginary 
part is plotted against the real part (Cole-Cole plot), x"(x') ( |Cole and Cole J 
1941; [Jonscher, 1983 ), where a classic Debye-relaxator should yield a perfect 



semicircle, centered on the x'-axis at (xo + Xs')/2 and with radius (xo — Xs')/2- 
The apex of the semi-circle corresponds to loTc = 1. Non-zero a has the effect 
to depress the semi-circle such that the angles between the x'-axis and the 
tangents at u = and — oo are =l=(l — a)7r/2, respectively. Figure ^ shows 
the susceptibility data for t = 0.9 (a) and 1.0 nm at different temperatures 
T = 45, 50, 55 and 60 K. The above derived expressions for the real and 
imaginary parts (Eq. |] and ^ can be expressed in the form ( [Hagiwara, 1998|) 



Xo - Xs 



2tan[(l - a)7r/2] \ 



+ \ ix' - Xs){Xo - X') + 



(Xo - XsY 



4tan2[(l - a)n/2]' 
(5) 

The fit yields similar results for a{T) compared to those from the fit to the 
imaginary part x"{f) (Fig- i|) solid versus open diamonds). It should be 
noticed that xs = in all cases, i.e. no measurable response is expected 
at frequencies above single particle flip frequencies. This corroborates the 
model of ferromagnetic order within each superparamagnetic particle. 



CONCLUSION 

The dynamical susceptibility of the SSG system [Co8oFe2o(t)/Al203(3nm)]io 
(t = 0.9 and 1.0 nm) was studied under the influence of a bias field and in 
view of its polydispersivity within the framework of a Cole-Cole description. 
Cole-Cole fits yield reasonable values for the characteristic relaxation time 
Tc of the system and for its polydispersivity exponent a. The relaxation 
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time can be well described by a critical power-law dependence. One should 
note that by this kind of extraction of Tc one has access to an extremely 
wide timescale of eight or ten orders of magnitude. Reasonably large values, 
a. «i 0.8, are obtained, which are typical of SG systems. The Cole-Cole plots 
of the susceptibility data confirm the SG characteristic, i.e. one observes a 
strongly flattened semi-circle. 
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Figure 1: x' and x" vs T ior t — 0.9 nm (a) and 1.0 nm (b) measured at 
constant frequency / = 0.1 (curves 1 and 1') and 1 Hz (curves 2 and 2') 

with ac field amplitude ^^Hac = 0.05 mT and dc field ^qH = mT (curves 
1 and 2) and ac field amplitude HoHac = 0.4 mT and dc field HqH = 0.6 mT 
(curves 1' and 2'), respectively. 
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Figure 2: Schematic drawing of the measurement conditions relevant for 
the data as numbered in Figure 1 (see text). The sohd curve shows M{H) 
without hysteresis. 
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Figure 3: x' {^) ^"^^ x" (b) vs / at different temperatures T = 45, 50, 55 
and 60 K for the t = 0.9 nm sample in zero-field and fioHac = 0.05 mT. Tlie 
lines in (a) are guides to the eyes and in (b) best fits according to Eq. ^. 
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Figure 4: x' {^) ^"^^ x" (b) vs / at different temperatures T = 45, 50, 55 
and 60 K for the t = 1.0 nm sample in zero-field and fioHac = 0.05 mT. The 
lines in (a) are guides to the eyes and in (b) best fits according to Eq. |. 
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Figure 5: Results from Cole-Cole fits to the data shown in Fig. ||, H and | 
for t = 0.9 (a) and 1.0 nm (b). The characteristic relaxation times Tc (open 
circles) are best fitted to a critical power law (solid and broken lines; see 
text). The polydispersivity exponent a vs T as obtained from fits to Eq. | 
(open diamonds) and to Eq. ^ (solid diamonds), respectively, are connected 
by eye-guiding lines. Encircled data points are results from non-converging 
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Figure 6: Cole-Cole plots of x" vs x' for t = 0.9 (a) and 1.0 nm (b) at 
different temperatures and frequencies as indicated. The solid lines are best 
fits according to Eq. ^ 



17 



